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CN Three-dimensional conformal field theories (CFTs) of deconfined gauge fields coupled to 

^ gapless flavors of fermionic and bosonic matter describe quantum critical points of condensed 

^ matter systems in two spatial dimensions. An important characteristic of these CFTs is the 

> finite part of the entanglement entropy across a circle. The negative of this quantity is equal 

^ to the finite part of the free energy of the Euclidean CFT on the three-sphere, and it has 

^ been proposed to satisfy the so called F-theorem, which states that it decreases under RG 

flow and is stationary at RG flxed points. We calculate the three-sphere free energy of non- 
supersymmetric gauge theory with a large number A^^^ of bosonic and/or fermionic flavors to 
the flrst subleading order in l/Np. We also calculate the exact free energies of the analogous 
chiral and non-chiral M = 2 supersymmetric theories using localization, and find agreement 
with the l/Np expansion. We analyze some RG fiows of supersymmetric theories, providing 
further evidence for the F-theorem. 
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1 Introduction 



Many interesting quantum critical points of condensed matter systems in two spatial dimen- 
sions [l 14 are described by conformal field theories (CFT) in three spacetime dimensions 
where massless fermionic and/or bosonic matter interacts with deconfined gauge fields. These 
include critical points found in insulating antiferromagnets and d-wave superconductors and 
between quantum Hall states. Such CFTs can be naturally analyzed by an expansion in 
1/Np, where Np is the number of 'flavors' of matter. This large Np limit is taken at fixed 
Nc, where is a measure of the size of the gauge group e.g. the non-abelian gauge group 
U (Nc). Classic examples of such CFTs include three-dimensional U (1) gauge theory coupled 
to a large number of massless charged scalars 15 or Dirac fermions [l6, 17 . These theories 
are conformal to all orders in the 1/Np expansion, and they are widely believed to be confor- 
mal for Np > Ncrit, where A'crit is a conjectured critical number of flavors dependent on the 



choice of the gauge group 16 17 . The 3-dimensional CFTs may also contain Chern-Simons 
terms whose coefficients k may be taken to be large. 

An important characteristic of a 3-dimensional CFT is the ground state entanglement 
entropy across a circle of radius R. Its general structure is 



S 



R 



a- 



where e is the short distance cut-off. As estabhshed in [T^ (see also 19,20]) the sub leading 
i?-independent term is related to the regulated Euclidean path integral Z of the CFT on 
the three-dimensional sphere S^: F = — log|Z|. The quantity F has been conjectured to 
decrease along any RG flow (l8|[2T|j23] p] This conjecture was inspired by the c-theorem in 



two spacetime dimensions 25 and the a-theorem in four spacetime dimensions 26 28 



In any 3-dimensional field theory with M >2 supersymmetry, the 5*^ free energy F may 



be calculated using the method of localization 29 32j. It has also been calculated in some 
simple non-supersymmetric CFTs, such as free field theories [T9l|20l|23l|33] and the Wilson- 



Fisher fixed point of the 0{N) model for large N 23 , which has been conjectured 34 to 



be dual to Vasiliev's higher-spin gauge theory in AdS^ (35). In this paper we present the 
calculation of F in certain 3-d gauge theories coupled to a large number of massless fiavors, 
to the first sub leading order in 1/Np. We will find that this subleading term is of order 
logA^i.. 

The CFTs we study have the following general structure. The matter sector has Dirac 



^After the original version of this paper appeared, a proof of the F-theorem was presented in 
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1 . . . Nf, and complex scalars, Za, a = 1 . . . Ni,. We will always take the 



fermions ipa, ct 

large Nj limit with N^/Nf fixed, and use the symbol Np to refer generically to either Nf or 
Nf,. These matter fields are coupled to each other and a gauge field by a Lagrangian of 
the form 



a=l a=l ^ ^ 



1.2) 



where = d^ — iA^ is the gauge covariant derivative, and the ellipses represent additional 
possible contact-couplings between the fermions and bosons, such as a Yukawa coupling. 
The scalar "mass" s generally has to be tuned to reach the quantum critical point at the 
renormalization group (RG) fixed point, which is described by a three dimensional CFT; 
however, this is the only relevant perturbation at the CFT fixed point, and so only a single 
parameter has to be tuned to access the fixed point. In some cases, such scalar mass terms 
are forbidden, and then the CFT describes a quantum critical phase. All other couplings, 
such as u and the Yukawa coupling, reach values associated with the RG fixed point, and so 
their values are immaterial for the universal properties of interest in the present paper. 

The gauge sector of the CFT has a traditional Maxwell term, along with a possible 
Chern-Simons term 



1 ^ „9 ik ^ 
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1.3) 



The gauge coupling has dimension of mass in three spacetime dimensions. It fiows to an 
RG fixed point value, and so its value is also immaterial; indeed, we can safely take the limit 
— oo at the outset. However, our results will depend upon the value of the Chern-Simons 
coupling /c, which is RG invariant. We will typically take the large Np limit with k/Np fixed 
at fixed and in most of this paper we set Nc = 1 for simplicity. (This is to be contrasted 
with the 't Hooft type limit of large Nc where k/Nc is held fixed; see, for example, recent 



work 36 38].) One of our principal results, estabhshed in section [3| is that for the U{1) 
gauge theory with Chern-Simons level fc, coupled to Nf massless Dirac fermions and A^j, 
massless complex scalars of charge 1 as in (1.2) with s = u = 0, 



^(Ar, + Ar,) + ^(Ar,-Ar,) + llog 




+ 



;i-4) 



This formula shows that the entanglement entropy is not simply the sum of the topological 



2 



contribution — | log k and the contribution of the gapless bulk modes, unlike in the models of 
For CFTs with interacting scalars relevant for condensed matter applications, we have to 



39 



consider the m — )■ oo limit, and this yields a correction of order unity, with F — )■ F— C(3) / (Svr^ 



at this order 123]. All higher order corrections to (1.4) are expected to be suppressed by 



integer powers of 1/Np, whose coefficients do not contain any factors of logNp. 

In section |4] we will examine similar Af = 2 supersymmetric CFTs on using the local- 
ization approach. We consider theories with chiral and non-chiral flavorings. The partition 
function Z on is given by a finite-dimensional integral, which has to be locally minimized 
with respect to the scaling dimensions of the matter fields [3ll|. For a theory with charged 
superfields we develop expansions for the scaling dimensions and for the entanglement 
entropy. As for the non-supersymmetric case, the subleading term in F is of order logA^. 
The coefficient of this term computed via localization agrees with the direct perturbative 



calculation (1.4) 



In the supersymmetric case it is possible to develop the 1/A^ expansions to a rather 
high order, and we compare them with precise numerical results. This comparison yields 
an unprecedented test of the validity and accuracy of the 1/A^ expansion. At least for 
supersymmetric CFTs, we find the 1/A^ expansion is accurate down to rather small values 
of A^. We also note a recent numerical study f40], which found reasonable accuracy in the 
1/Nf) expansion for a non-supersymmetric CFT. 



2 Mapping to and large Np expansion 

Let us start by examining the case of a f/ (1) gauge field. After sending — )■ oo, the combined 



Lagrangian Cm + Ca obtained from (1.2) and (1.3) contains two relevant couplings s and u 



and we should first understand to what values we need to tune them in order to describe an 
RG fixed point. Let's ignore for the moment the fermions and the gauge field and focus on 
the complex scalar fields. The path integral on a space with arbitrary metric is 



With the help of an extra field A, this path integral can be equivalently written as 



(2.1) 



Z = C DzaDX exp 



(2.2) 



3 



where the normahzation factor C defined through C f DXexp [— f d?r ^ (2^-^^)] — 1 
introduced so that the value of the path integral stays unchanged. 

In fiat three-dimensional space, we can tune s = m = and describe a non-interacting 
CFT of A'^^ complex scalars. If instead we tune s = and send u — ?■ 00, the path integrals 



(2.1 ) and (2.2) describe the interacting fixed point that we will primarily be interested in in 
this paper. We can also send both s and u to infinity, in which case the path integrals above 
describe the empty field theory. Using conformal symmetry, we can map each of these fixed 
points to by simply mapping all the correlators in the theory. Indeed, since the metric 
on is equal to that on up to a conformal transformation. 



53 



(1 



\r\ j 



(2.3) 



the mapping of correlators to S*^ is achieved by replacing 



2ff - f 



r — r ^ 



(2.4) 



in all the fiat-space expressions]^ While the theory on defined this way certainly has 
the correlators of a CFT, it may be a priori not clear which action, and in particular which 
values of s and u, one should choose in order to reproduce these correlators. 

In order to study the free theory on 5''^ one should tune s = 3/4 and m = 0. This result 
holds to all orders in and one can understand it as follows. The two-point connected 
correlator of Za on M? is 



{Za{r)zb{r')) 



Jab 



free 



(2.5) 



because it is the unique solution to the equation of motion following from (2.1 ) with a delta- 
function source, — V^a (^a('")^fe('"'))frce — SabS^^\r — r'). Using the mapping (2.4) we infer 
that the corresponding two-point correlator on 5*^ should be 



(Za(r)2b(r ))fee 



Sab (l + |rf)'/'(l + 




2^1/2 


Stt 







(2.6) 



^This replacement certainly works for correlators of scalar operators. In the case of vector operators it 



is still true that one can use (2.4 1 provided that the correlators are expressed in a frame basis, as in the 
following section. 
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An explicit computation shows that (— V^s + 3/4) {za{r)zb{r'))f^l^ = SabS^^\f — r')/ ^Jg{r), 



which is indeed the equation of motion that would follow from (2.1) with s = 3/4 and 
u = 0. This result was of course to be expected because a mass squared given by s = 3/4 
corresponds to a conformally coupled scalar. 

A more subtle issue is how to map to the interacting fixed point, which in fiat space 
had s = and ^ = 00. As explained for example in j4l], the generating functional of 
connected correlation functions of the singlet operator \za\'^ in the theory with u = 00 equals 
the Legendre transform of the corresponding generating functional in the theory with m = 0, 
to leading order in a large A*";, expansion. This result holds on any manifold, and in particular 
on both M'^ and S^, and it assumes the other couplings in the theory are held fixed. If we 
set s = on M'^ and s = 3/4 on S^, the Legendre transform assures us, for example, that to 
leading order in A^;, the two-point correlators in the theory with u = 00 are 



(l^.('-)Pk.('-')RS«c., = 7Z^. 



2\2 /, , |^|2\2 (2.7) 



16 \r — r'\ 

with the same normalization constant c, which is consistent with the conformal mapping 



of correlators realized through eq. (2.4). While in the free theory Za is a free field and the 
operator \za\^ therefore has dimension 1, in the interacting theory | Z(i I is cL dimension 2 
operator. To study the interacting fixed point on we therefore should set s = 3/4 + 



0{1/Nb) and take m 00 in (2.2). 



Reintroducing the fermionic and gauge fields, we can write down the action as 

S = j dh^ ^a^'^D^ij^ + \D^Za\'' + is- zA) \Za\^ + + J A A dA . (2.8) 

This action is of course invariant under gauge transformations, and therefore a correct defi- 
nition of the path integral requires gauge fixing: 

where Vol(G) is the volume of the group of gauge transformations, and X denotes generically 
all fields besides the gauge field. One justification for this normalization of the path integral 
is that for a pure Chern-Simons gauge theory on it yields the expected answer |42| 
Z = l/\/k, as will emerge from our computations below. Because the first cohomology of 



5 



is trivial, we can write uniquely any gauge field configuration A as A = B + d(f), where 
d*B = and (p is defined only up to constant shifts. One should think of B as the gauge-fixed 
version of A and of d(f) as the possible gauge transformations of A. Since the action S[A, X] 
is gauge- invariant, it is independent of and only depends on B: S[A,X] = S[B,X]. 

We claim that 



DA = DB D{d(j)) = DB D'cP^det ' (-V2) , (2.10) 

where D'(j) means that we're not integrating over configurations with = constant, and det' 
denotes the determinant with the zero modes removed from the spectrum. To understand 
this relation, first note that the space Qp{S^) of p-forms on is a metric space with the 
distance function V(u,uj + 5co) = (/ \5ujf) ^ . Then V{d(f),d(f) + d5(j)) = (/ ^ = 

(J S(f) (— V^) 50)"^^^ after integration by parts, and also T'(0, (p + 5(f)) = (J |50|^) • In other 
words, for each component of in a basis of eigenfunctions of the Laplacian, the distance 
between d(j) and d(j) + d6(f) is larger than the distance between (p and + 50 by a factor of the 



square root of the eigenvalue with respect to — V^. Eq. (2.10) follows as a straightforward 
change of variables. 

The gauge transformations are maps from into the Lie algebra of the gauge group. 
The volume of the group of gauge transformations Vol(G) can be expressed as 



Vo\{G) = Vo\{H) j D'(j), (2.11) 



where H is the group of constant gauge transformations, and J D'cf) is an integral over 
the non-constant gauge transformations with the measure given by the metric function V 
introduced in the previous paragraph. In the case of a compact f/(l) with Vol(f/(l)) = 27r, 
a constant gauge transformation = c has c G [0, 27r). Therefore 




Combining (2.9)-(2.12) we obtain 



27r/Vop3) J 



(2.13) 



In this paper we will use the partition function in eq. (2.13) to compute F = — log \Z\ in the 



6 



limit where Nf, Nf,, and k are taken to be large and of the same order. 

To leading order in the number of flavors we can ignore the gauge field and the Lagrange 
multiplier field A. Setting s = 3/4 as discussed above, we can write down the resulting path 
integral as 



Dipa Dza exp 



(2.14) 



In this approximation we have a theory of free Nf Dirac fermions and A'^^ complex scalars 



with the free energy 23 



i^o = ^(iV/ + iV,) + ^(iV,-iV,) 



(2.15) 



To find the corrections to Fq we write (2.13) approximately as 
^^,-i^oCv/det'(-V^) 



27rv/Vol(S3 



(2.16) 



with 



SA>^]= I dh-^)^\{rf-\ I d'ryW) f rfVv^A(r)A(r')(k.(r)n^,(r')0,l 



where 



S3 

free ' 



(2.17) 



J/^(r) = i)a{r)-f''ipa{r) + iza{r)d^ Za{r) - iZa{r)d^'Za{r) . 



(2.18) 



In writing the effective action (2.17) we used (|-2a('")P)iee — 0; which follows from the fact 



that the free theory < \2.1A\ is a CFT. 
Defining 



5Fx 
SFa 



log 
log 



C j DAe'^^«[^l 
27rx/Vol(S3) J 



(2.19) 



7 



we can then write F = Fo + 5Fx + 5 Fa + o{N°). The quantity 5 Fx was computed in 



23 



6Fx = -^^. (2.20) 



We devote the next section of this paper to calculating SFa- 



3 Gauge field contribution to the free energy 



3.1 Performing the Gaussian integrals 

Let's denote by K^'^ the integration kernel appearing in S^, namely 

ik 6^{r — r') 1 



-{nr)r{r')) 



free ' 



2tt 



9[r) 



■g{r') " 



(3.1) 



As discussed above, when one writes A = dcf) + B the action should be independent of 0, 
so pure gauge configurations Ay{r') = dl(f){r') are exact zero modes of the kernel K'^'^{r, r'). 
Since we should integrate over the gauge-fixed field B only, the Gaussian integration of the 
effective theory yields 1 / ^/det^{IO^^^J{2^ . A gain, the prime means that we remove 

the zero modes from the spectrum when we calculate the determinant. 

Reinstating the radius R of the sphere measured in units of some fixed UV cutoff, the 
discussion in the previous two paragraphs can be summarized as 



SFa 



tr' log 



tr' log 

2 ^ 



log (2n^/R^Vo\{S^ 



(3.2) 



where all the operators are defined on an of unit radius. Out of the first two terms in 
this expression, the second one is the easier one to calculate (see also |4^). The spectrum of 
the Laplacian on a unit-radius consists of eigenvalues n{n + 2) with multiplicity {n + 1)^ 
for any n > 0. One first rearranges the terms in the sum as 

1/1 f ^^2 1 n{n + 2) 2 n log(2/i?) , , 

2 tr' log f - — j = - 5^(n + 1)2 log = J^i^' + 1) log ^ - . (3.3) 

^ ^ 71=1 n=l 



Then, using zeta-function regularization one writes 



tr' log 



i?2 J 



log(2/i?) d 



ds 



E 



+ 1 



{ (n/RY 



C(3) log(7ri?2) 

id + • (3-4) 



8 



Combining this expression with eq. (3.2) and using Vol(S'^) = 27r^, we obtain 



8Fa = ^ tr' log 



C(3) logjSn'R) 
4772 ^ 2 



(3.5) 



The only remaining task is the computation of the first term in this equation that we perform 
in the next subsection by explicit diagonalization of K^^". 

3.2 Diagonalizing the kernel K^^^{r,r') 

Ultimately we would like to diagonalize the kernel K'^'^ on S^. However, as a warm up it is 
instructive to consider the same diagonalization problem in flat space first. 

3.2.1 Warm-up: Diagonalization on 

The first step is to calculate the two-point function of current {J^ (r) J'^ {0))^l^, where we use 
the superscript to emphasize that we are in flat space. If we normalize the two-point 
functions of Za and ipa to be 



ab 



Jab 



(27r)3|p|2 47r|r| 



(3.6) 



|3 ' 



(27r)3 47r \r\ 

then the two-point function of the current may be straightforwardly calculated to be 



Nf + Nb\r\ SI"" -2r''r'' 



free 



(3.7) 



It is simple to check that eq. (3.7) is of the right form. This correlator is fixed up to an 



overall constant by the requirements that it should be homogeneous of degree —4 in r [J^ is 
a dimension 2 operator) and that away from r = it should satisfy the conservation equation 
5m(^^(^)^'(0))m3 = 0. Using 



-vr' \p\ , 



d-^r 



Jp-r ^2 



12 



(3.^ 



and introducing the Fourier space representation of the kernel K^'^ via 



K^,u{r, r') 



(27r)= 



(3.9) 



9 



one obtains 
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KM = \P\ (^.. - Ip') + ^'^^'P^ ■ (3-10) 



For fixed p, the eigenvalues of this matrix are 



0, MfMi.i, (3.11) 



2 V 8 TT 

The eigenvector associated with the zero eigenvalue is as expected ipi,e^^'^'' , corresponding to 
a gauge configuration A^, = d^cf). We will now see that on S*^, while the diagonalization of 
K^'^ is significantly more complicated, the answer is equally simple: the magnitude of the 



momentum p appearing in (3.11) should be replaced by a positive integer label n. 



3.2.2 Diagonalization on 

When we work with vector fields on it is convenient to introduce the dreibein 

2 



e r 



1 I I |2 

1 + \r\ 



dr' (3.12) 



and work only with frame indices. For example, 

{Ar)J={r'))L = e',{r)ei{r'){r{r)r{r'))ll . (3.13) 

The frame indices i and j are raised and lowered with the flat metric, so there is no distinction 
between lower and upper frame indices in Euclidean signature. 



Using the transformation of correlators under Weyl rescalings in eq. (2.4), one can im- 
mediately write down the current two-point function on S^: 

iJ'ML-'^^^^'^^^^. (3.14) 

As in fiat space, the form of this correlator is fixed by the requirement that away from r = 
we must have V i{J\r)J^ {Q)) = 0. 

To understand the diagonalization of K"^^ on we need to know that the space of square- 
integrable one-forms on S^, being a vector space acted on by the S'0(4) = SU{2)l x SU{2)r 
isometry group, decomposes into irreducible representations of 5*0(4) as follows. Any one- 



10 



form bj can be written as a sum of a closed one-form and a co-closed one-form. The closed 
one-forms on 5''^ are of course cohomologous to zero, so they're also exact. A basis for 
them therefore consists of the gradients of the usual spherical harmonics. Like the spherical 
harmonics, they transform in irreps with ji = j^. On the other hand, the co-closed one-forms 
transform in irreps with = j^^l. So an arbitrary one- form can be written as 

n,£,m 

where we denoted by S"^*" the closed component transforming in the irrep with ji = jn = 
{n — l)/2 and by V^^™ and V^^™ the co-closed components transforming in the irreps with 
Jr = Jl + ^ = n/2 and ji = Jr + ^ = n/2, respectively. All the harmonics appearing in 



( |3J[5| have n > 2. For §f there are n states in each irrep indexed by the integers i and 
m satisfying < i < n and —i<m<i. For the other two classes of vector harmonics, we 
have the same bounds on m but now < i < n, giving a total dimension of — 1 for each 
irrep. 

The 5*0(4) generators commute with the kernel Kij, so the eigenvectors of this kernel 
can be taken to be S"^™, V^^™, and V^^T*- The spectral decomposition of Kij is therefore 



K, 



n,£,m 

where Sn, v^, and are the corresponding eigenvalues. These eigenvalues are independent 
of i and m because for fixed n one can change i and m by acting with the 30(4) generators, 
which commute with Kij. The degeneracy of s„ is and that of and is — 1, with 
n>2. 

Given Kij one can find its eigenvalues by taking inner products with the eigenvectors. 
Using rotational invariance, one can actually set r' = after summing over £ and m. For 
example. 



_ Vol(^3) 

5 77, 



J d\Sf"'{ryKij{r,0)Sf"'{0), (3.17) 



^2 z — ^ z — ^ (3 

£=0 m=-£ 



where the n in the denominator is the dimension of the 5*0(4) irrep to which belong 



11 



We notice that only the harmonics with i = 1 contribute, so 



Vo\{S^) V 



V [ rfV§f™(r)*ir,,(r,0)§f™(0) 



(3.18) 



m=— 1 



with similar expressions for and v^, the only difference being that n"^ should be replaced 
by — 1. 

Using explicit formulae for the harmonics (see Appendix |A]) , one obtains 

Nf + Nh 



647rn(n2 — 1) 



/ dx csc^ — sin x [—2n cos(nx) sin x+ — n'^ + cos xin'^ + 1)) sin(nx)] 
Jo 2 

(3.19) 



and 



r p A^f + Nh r-x r / \ / 2 2 \ 1 

vj = — — -^-^ — — I ax CSC — sin X p sin X cos(nx) + — n cos x — Ij sinnxj ± 



QA'Kn{rfi — 1) Jq 



27r 
(3.20) 



The integration variable x appearing here is related to r through |r| = cot(x/2). 



We expect s„ = because of gauge invariance. Both (3.20) and (3.19) are divergent at 



X = 0, and need to be regulated. A way of regulating them that doesn't preserve gauge 
invariance is to replace csc^ | by esc" |, compute these integrals for values of a for which 
they are convergent, and then formally set a = 6. Another way would be to assume S2 = 0, 
and calculate Sn — S2 and Vn — S2, which are now convergent integrals. Both of these ways 



of regulating (3.19) and (3.20) give 



Sn = 0: 



^ _ n{Nf + Nb) ikn 
~ 16 ^ 27 



^ _ n{Nf + Nb) ikn 



(3.21) 



16 2n 

Note the similarity between these expressions and the corresponding flat-space ones in 
eq. (KU). 
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3.3 Contribution to the free energy 

We can now evaluate the first term in (3.5): 



- tr' loe f 

2 ^ V 2'kR 



^(n2-l)log 



n=2 
1 



2tiR 



log 



+ 



m_\ogR 

47r2 2 



(3.22) 



where the second line was obtained with the help of zeta-function regularization. Combining 



this expression with (3.5) yields 



5Fa 



log 




(3.23) 



Note that all of the logi? terms cancel in the final answer, as they should since we are 
describing a conformal fixed point, for which the path integral should be independent of R. 
Another check of this result is that when Nf = N), = we recover the standard result for 
the free energy of U (1) CS theory on 42 , 6Fa = ^ log k. 



As an aside, we note that if we included the Maxwell term in (1.3), eq. (3.22) would be 
modified to 



- tr' log (^^^ = y^(n^ - 1) log ^— (v^ + 
2 ^ \2ttR y ^ 27iR \ " 



n 



(3.24) 



with still defined as in (3.21). Of course, flows to infinity in the IR, so as long 
as we have a non-zero CS level or non-zero numbers of flavors one can safely ignore the 



contribution from the Maxwell term in (|3.24|). If however one studies pure Maxwell theory 
with k = Nf 



so that v^ = Om ( |3^ , the free energy become^ 

\og{e'R) , C(3) 



Maxwell 



+ 



4772 



(3.25) 



The logarithmic dependence on R is consistent with the fact that the free Maxwell theory 
is not conformal in three spacetime dimensions. i^Maxweii decreases monotonically from the 
UV (small R) to the IR (large R). 

^We thank D. JafFeris and Z. Komargodski for very useful discussions of the free Maxwell field on S^. 
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3.4 Generalization to U{Nc) theory 



Eq. (3.5) generalizes straightforwardly to the case of U{Nc) gauge theory with Nf Dirac 
fermions and A^^ complex bosons transforming in the fundamental representation of the 
gauge group. At large k the contribution of the Tr term in the Chern-Simons Lagrangian 
(1.3) to the partition function is suppressed by 1/ y/k, and the quadratic term proportional 



to Tt A A dA is the same as that of A*"^ U{1) gauge fields with Chern-Simons coupling k. 
There are Nf Dirac fermions and A^b complex bosons charged under each of these U (1) gauge 



fields. One then just has to multiply the U{1) answer (3.23) by a factor of A"^: 



5Fa 



iV2 




+ log 



Vol([/(A/'e)) 

Vol(f/(l))^c 



(3.26) 



The second term in this expression comes from the different gauge fixing of the U {N^) gauge 
theory compared to a theory of A"^ f/(l) gauge fields. As explained in section [2| the gauge 
fixing procedure involves dividing the partition function by the volume of the gauge group, so 
the partition function for the U{Nc) theory has a prefactor of 1/ Vol(?7(A"c)) while the f/(l)^= 
theory obtained by multiplying (3.23) by A"^ would have a prefactor of 1/ Vol(f/(l))^^. We 
have (see for example [43]) 



V0l(t/(Ar,)) (27,)-^c(^c-l)/2 



Vol(?7(l))^c l!.2!---(Ar,-l)! 



(3.27) 



Thus, for U{Nc) gauge theory with Nf fundamental fermions and A^;, fundamental bosons 
we have 




Nf + NA^ , rk^ ' 

vr 



--N,{N, - 1) log(27r) - log (1! ■ 2! ■ ■ ■ (AT, - 1)!) + . . . 



(3.28) 



with corrections expected to vanish in the limit of large Np. In writing (3.28) we kept N^. of 
order one while scaling Nf, N^, and k to infinity with their ratios fixed. Generalizing (3.28) 
to different gauge groups proceeds in a similar way. 
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4 SUSY gauge theory with flavors 



In this section we compute the free energy of U (1) Chern-Simons matter theories with M >2 
supersymmetry coupled to a large number of flavors. These computations allow us to check 



the first sub-leading correction to the non-SUSY result in the equation (3.23) in a different 
way. The computations in this section have as starting point the results of refs. 29 ,31 , 
which used the technique of supersymmetric localization to rewrite the partition function 
of theories with M >2 SUSY as finite-dimensional integrals. Our computations also involve 
finding the scaling dimensions of the gauge-invariant operators. 

4.1 A/" = 4 theory 

As a warmup to the M = 2 calculations, consider the A/" = 4 parity-preserving supersym- 
metric f/(l) theory consisting of A/" = 4 hypermultiplets coupled to an A/" = 4 vector 
multiplet. In A" = 2 notation, the A/" = 4 vector multiplet consists of an A/" = 2 vector and 
a neutral chiral superfield $ of dimension 1. The A/" = 4 supersymmetry does not allow a 
Chern-Simons term. The hypermultiplets can be rewritten as pairs of oppositely charged 
chiral-multiplets Qa of U{1) charge +1 and Qa with U{1) charge —1. The Af = 4 SUSY 
requires a superpotential interaction W ~ Qa^Qa- The superpotential has R-charge equal 
to 2. Then the SU{2) subgroup of the S0{4)r R-symmetry, under which Qa and Qa trans- 
form as a doublet, fixes the R-charge of the matter chiral multiplets to have the canonical 
free-field value: Aq = Aq = 1/2. The partition function is then given by 



29 



d\ 2-^r (f ) 



2^ cosh^(7rA) v^r (^) 
Expanding this at large we find 



(4.1) 



F=-logZ = A'log2 + llog(^)-i^ + ^ + .... (4,2) 

This large expansion is asymptotic, but it provides a very good approximation of the 
exact answer (4.1) down to A^ = 1 — see Figure [!} Including more terms in (4.2) makes the 



approximation worse at A^ = 1. 

With A^ pairs of hyper-multiplets we have a total of 2A^ physical complex bosons and 2A^ 



Dirac fermions. We then see perfect agreement of the first two terms in (4.2) with eqs. (2.15) 
and (|333|. 



15 



10 



N 



Figure 1: The exact free energy of the A/" = 4 theory obtained from eq. (4.1) (sohd orange 



and the analytical approximation (4.2) (dashed black) 



4.2 A/" = 3 theory 

Let us add the Chern-Simons term for the M = 2 abelian vector multiplet; it breaks A/" = 4 
down to A/" = 3 supersymmetry. The field content is the same as that of an A/" = 4 vector 
multiplet and hypermultiplets, namely an A/" = 2 vector, a neutral chiral $, and pairs 
of chiral multiplets Qa and Qa charged under the M = 2 vector. The superpotential required 
by A/" = 3 SUSY is 



W 



k 

Air 



(4.3) 



After integrating out the massive field $, the superpotential can be rewritten as 44 



W 



27r,^ 



k 



{QaQaY 



(4.4) 



The conformal dimensions of Qa and Qa are still equal to 1/2, as is required by the marginal- 
ity of W and by the Z2 symmetry under which Qa and Qa are interchanged and all the fields 



in the vector multiplet change sign. The partition function is 29 



1 

2^ 



dX- 



cosh^(7rA) 



(4.5) 
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While this expression cannot be evaluated analytically, one can evaluate it using a saddle 
point approximation in the limit where both k and N are taken to be large. Let us define 
K = 2k/{N7[) and take to infinity while keeping k fixed. The saddle point is at A = 0, 
and in order to obtain a systematic expansion, one should write 



1 

2^ 



1 + 



12 



45 



+ 



10080 



+ 



(4.6) 



where the parenthesis contains the small A expansion of the function e^^'^'/^ cosh'^ (nX). 



Order by order in this expansion one can perform the integrals in (4.6) analytically. The 
result is 



2^ y A^7r(l - zk) 



1 



1 



68 + 35iV 

4A^(1 -m)2 ~ 3iV2(l ^ 96iV3(l - zk)^ 



+ 



(4.7) 



Calculating F = — log \ Z\ and expanding in N, we obtain 



F = iVlog2+-log( 



Nn 



4/«2(/c2_l) 



+ 0{N- 



(4. 



We note that in order to calculate the 0{N ") term in F we need the expansion in (4.6) to 



be up to order 0(A^"). The expression ( |4.8 ) is also in agreement with eqs. (2.15) and (3.23) 
given that the J\f = 3 theory has Nh = Nj = 2N. 

Let us extend our discussion to the non-abelian theory with gauge group U{Nc). The field 
content now consists of an = 4 vector multiplet in the adjoint representation of the gauge 
group and pairs of chiral multiplets Qa and Qa, in the fundamental and anti-fundamental 
representations of the gauge group, respectively. After localization the partition function for 
this theory is given by [29] 



2N4Nc-l) 



Y[dX,] I JJsinh2[7r(A, - A^)] ) exp ( tnkj^^^i ] Y[cosh-^ (nXi) . 

(4.9) 



,i=l 



i=l 



i=l 



In the limit where N^/N <^ 1, the integral has a saddle point at Aj = O [(A'^^/A^)^''^] . 
Through next to leading order the partition function of the non-abelian theory reduces to 



(27r)^-(^=-i) 



NAN- 



i=l 
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where k is defined as in the abehan theory and Aj = v N\i. We have rescaled the integration 



variables so that the remaining integrals in eq. (4.10) produce numbers independent of 



Taking the log of eq. (4.10) we then see immediately that 



Ar2 

F = N^N \og2 + ^log(A^) + 0{N^) . 



(4.11) 



Given that the non-abelian theory has A^^ = N j = 2N, the equation above is in agreement 



with eq. (3.28). 



4.3 Non-chiral J\f = 2 theory 

Moving up one notch in complexity, we now consider the M = 2 Chern-Simons theory coupled 



to the chiral fields Qa and Qa introduced above, this time without the superpotential (4.4) 



The absence of the superpotential leaves the R-charges of Qa and Qa a priori unrestricted. It 
was proposed in [31] that one way of finding the correct IR R-charges in an A/" = 2 theory is by 
calculating the partition function on for any choice of trial R-charges consistent with the 
marginality of the superpotential and then extremizing over all such R-charge assignments. 
The R-charges of Qa and Qa can be taken to be equal to some common value A because of 
the following symmetries: the action is invariant under two U{N) symmetries under which 
the Qa and Qa transform as fundamental vectors, as well as under a charge conjugation 
symmetry that flips the sign of all the fields in the vector multiplet and at the same time 
interchanges Qa and Qa- 



As a function of A, the partition function is 31 



dXe' 



TfcA2 Af(^(l-A+iA)+£(l-A-iA)) 



(4.12) 



where the function i{z) is given by 



■ log (l 



2TTiz 



) + 



TCZ 



+ ^Li2 {e'-") ] - 



12 



(4.13) 



This function can be found by solving the differential equation dz^{z) = —iTzcot{TTz) with 
the boundary condition i{0) = 0. It is a real function when z is real. 

We again take to infinity while keeping k = 2k /{Nn) fixed. In this limit one can use 



the saddle point approximation to calculate the partition function (4.12) as in the previous 
section. The exponent in (4.12) is an even function of A, so there is a saddle point at A = 0, 
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and we will assume this is the only relevant saddle. To leading order in we therefore have 



F(A) = -2N^{l - A) + 0{\ogN) . 



(4.14) 



This function is maximized when dF/dA = 27r(A — l)cot(7rA) = 0, which implies A = 
1/2 + 0{N~^). We will find that this anomalous dimension affects F only at order 1/A^, i.e. 
the first two leading orders in the large N expansion of F are the same for the M = 2 theory 
and the A/" = 3 theory studied in the previous section. 

One can develop a systematic expansion to study 1/A^ corrections in a similar way to 
what was done at the end of the previous section for the A/" = 3 theory. The fact that now 
A depends on N introduces an extra complication. We expand A as 



1 Ai As 



(4.15) 



and we rescale A = X/yN. One can then write 



2^ViV 



evr^A? + 24AiA2 + 

1 H h . . . 

UN 



(4.16) 



where the expansion in parenthesis is in powers of 1/A^ while holding A fixed. Term by term 
in this expansion, these integrals can be evaluated analytically. The free energy is 



1 fN-K 
F(A) = iVlog2+-log( — 



TT^A? 2Ai 1 - 



1 + ^2 4(1 + k2)2 J N 



1 



+ 



(4.17) 



Maximizing this expression with respect to Ai we obtain 

A.: 



7r2(l + ^2 



(4.18) 



For ^ ^ 1 this result agrees with section 6.3 of 31 . Repeating this procedure two 
more orders in F we find 



A 



1 



1 2 [7r2 - 12 + /€2(4 _ 27r2) + 7r2/c4] i 



2 tt^(1 + k'^)N 



7r4(l + K 



2\3 



iV2 



+ 0{N-'). 



(4.19) 



This series appears to be perfectly convergent. In fig. [2] we plot A(A^) for a few values of n 



using both the precise numerical result and the approximation (4.19). 
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Figure 2: The R-charge A plotted as a function of N for k = 0,4/7r,8/7r, with darker 
plots corresponding to larger k. The solid lines are calculated using the approximation in 



eq. (4.19). The circles are computed by numerically maximizing the free energy with respect 



to A. Note that the two computations match well even for small N. 



Using eqs. (4.17) and (4.18) we find that the free energy is 



+ 



+ 0{N-^). (4.20) 



Using Nf, = Nf = 2N, we see that this expression agrees with eqs. (2.15) and (3.23) that were 



derived directly from a large N expansion without the use of supersymmetric localization. 
Let us perturb the Af = 2 theory discussed above by the quartic superpotential 



W = g{QaQa) 



(4.21) 



Since, as can be seen from (4.18), the dimension of Qa and Qa is slightly smaller than 1/2, 
the perturbation (4.21) is a slightly relevant perturbation of the UV M = 2 theory. This 



theory should flow to an IR fixed point where the superpotential is exactly marginal, i.e. 
the IR R-charges of Qa and Qa are 1/2. The calculation of Fir is thus exactly the same as 
for the = 3 superconformal U{1) theory discussed in section 4.2. The infrared theory is 



conformal for any and for the special value = 27r/A; it is the A/" = 3 theory in eq. (4.4) 



Eqs. (4.8) and (4.20) imply that the change in free energy between the UV and IR fixed 
points is 



'UV 



'IR 



7r2 1 + ^2 2iV 



+ 0{N- 



(4.22) 
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which can be exphcitly seen to be positive, in agreement with the conjectured F-theorem [22]. 



Since the superpotential deformation (4.21) is only shghtly relevant, one may wonder 



how the result (4.22) compares with the perturbative computation performed in p3|. In [23] 
it was shown that if the Lagrangian is perturbed by a slightly relevant scalar operator of 
dimension 3 — e, then there is a perturbative IR fixed point and Fyy — Fin cc e^. If however 
the Lagrangian is perturbed by a pseudoscalar operator of dimension 3 — e, then there is no 
perturbative fixed point; it was seen in an example that if a fixed point exists then one might 



expect — F^ oc e. In our case, the superpotential deformation (4.21) translates into 



perturbations of the Lagrangian by both a scalar operator Oi and a pseudoscalar operator 
O2. Indeed, denoting 



Qo 



V2e^a + O^Fa , Qa = ^a + V^^^a + 0^ Fa 



(4.23) 



we have 



6C = g'Oi+g02, 

2 



01 



')a<Pa<Pb 



(4.24) 



The scaling dimensions of these operators are 



A(Oi 



6^ + 0(iV-2), A(02 



3 + 4 



Ai 



0{N' 



(4.25) 



so the pseudoscalar operator O2 is the more relevant one. One might expect the IR fixed 
point should be non-perturbative and that Fuv — -Fir oc — Ai/A^ times a function of order 
one. That the IR fixed point is non-perturbative can be seen after writing g = g/N so that 
g stays of order 1 as we take to infinity. The IR coupling giji = 2TT/k corresponds to 
gin = A/k, which is of order one in the large limit, meaning that the IR fixed point is non- 
perturbative. That FuY — -Fir oc —Ai/N times a function of order one can be immediately 



seen from eqs. (4.22) and (4.18). 



4.4 Chiral Af = 2 theory 

We now consider a natural generalization of the non-chiral M = 2 theory discussed in the 
previous section — the chiral M = 2 theory. This theory is given by A/" = 2 Chern-Simons 
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theory coupled to chiral fields Qa and anti-cliiral fields Qa with no superpotential. 
When N = N this theory reduces to the non-chiral theory discussed in the previous section. 
Without loss of generality, we now assume that N > N. Instead of dealing with N and 
it is convenient to define the following quantities, 

N = , /i = ^ , < /i < 1 . 4.26 

The R-charges of Qa and Qa, which we denote by A and A, are not gauge-invariant 
observables. Gauge invariant operators may be constructed from combinations of Qa, Qa, 
and the monopole operators Tm, which create m units of magnetic flux through 2-spheres 



surrounding their insertion points. The R-charge of Tm is given by 45 , 46 



R[Tm]=lim\+ni6, (4.27) 

where 7|m| is determined in terms of A and A, while 6 is so far arbitrary. In the F- 
maximization procedure one finds that in the space of S, A, and A there is exactly one 
fiat direction: F remains unchanged if we send simultaneously A— j-A + r, A— i-A — r, 
and S ^ S + kr for any r [22]. The R-charges of the gauge- invariant operators are of course 
independent of r. As long as A; 7^ 0, we can set 5 = as a gauge choice and work only with 
A and A, which are not necessarily equal when /i 7^ 0. 

As a function of the R-charges A and A, the partition function we need to consider is 
then 

^ink\2 ^Ne(l-A+i\)+Neil~A~iX) 23) 

We want to calculate the partition function in the limit where N goes to infinity and k = 
2k/{NTT) and /i are held fixed. In the large N limit we again find a saddle point at A = 0. 
The saddle point equation requires A = 1/2 + 0(l/iV) and A = 1/2 + 0{1/N). In order to 
study 1/N corrections, we expand the R-charges as 

1 Ai A2 ~ 1 Ai A2 I A 

A = - + ^ + ^+ A = - + ^ + ^+ (4.29) 

Using the methods developed in the previous sections, we can calculate the free energy 
perturbatively in the 1/N expansion and maximize the resulting expression term by term 
with respect to the Aj and Aj. Going through the procedure we find the following results 
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for the free energy and the R-charges: 
1 2(1 + 



A 



7r2(l + K2)iV 



+ 0{N- 



A 



1 2(1 -/x) 



2 7r2(l + K2)]V 



+ 0{N- 



F = Nlog2 + Uog(^VT 



1 



+ 



_4(l + /«2)2 n^l + K 



2\2 



(4.30) 



7r2 \^(1 + k2)2 3(1 + ^:2)3 



Using Nb = Nf = 2N, we see that the expression for F agrees with eqs. (2.15) and (3.23) 
that were derived without the use of supersymmetric locahzation. The combination A + A, 
which gives the R-charge of the the gauge invariant meson operators QaQt, is in agreement 



with the pertrubative calculations in 47,48 



We can perturb this theory by adding in the superpotentiaQ 



a.b 



(4.31) 



Since A + A < 1 in the UV M = 2 CFT, this superpotential deformation is relevant and 
causes an RG flow to the flxed point where the superpotential is exactly marginal. At the 
IR = 2 flxed point we have the constraint A + A = 1. To determine the free energy at 
the IR flxed point we simply have to repeat the F-maximization procedure above subject 
to this constraint. That in the UV one has to maximize F without any constraints while in 
the IR one has to maximize F under the constraint A + A = 1 means that the free energy 
of the IR flxed point is necessarily at most equal to the free energy of the UV flxed point. 
Indeed, we flnd that 



UV 



^IR 



2(1-^2 



1 



7r2(l + /t2)2jY 



0{N- 



(4.32) 



which is manifestly positive when /i2 < x When /i = 1 there are no Qa flelds, and so we are 
not allowed to add in the superpotential deformation. The R-charges at the IR flxed point 
are given by 



A 



4/x 



7r2(l + ^2) 



+ 0{N~ 



A = 1 - A. 



(4.33) 



^Changing the relative coefficients of the terms in (4.31 1 is an exactly marginal deformation [49j and does 
not change F. 
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5 Discussion 



In this paper we have studied certain 3-dimensional gauge theories coupled to a large number 
Np of massless charged fields. Such theories are conformal for a sufficiently large Np, and a 
good tool for studying them is the 1/Np expansion. In this paper we used such an expansion 
to study the disk entanglement entropy, which is related to the free energy F on the 3-sphere. 
For the U{Nc) gauge theory coupled to Nf massless Dirac fermions and Ni, massless scalars 



we found the first subleading term in the expansion, (3.28). We have also studied the M = 2 
supersymmetric abelian gauge theory coupled to positively charged chiral superfields Q 
and N negatively charged chiral superfields Q. In this case, F can be calculated numerically 
for any using the methods of localization. We compared these numerical results with their 
expansion and found excellent agreement down to small A^. 
An important question concerning such CFTs is whether there is a breakdown of confor- 
mal invariance for sufficiently small Np. In the M = 2 supersymmetric U{1) gauge theory, 
even for a single non-chiral flavor the theory is conformal and unitary. This is indicated by 



the mirror symmetry arguments 50 and conflrmed by explicit calculation of the localized 
path integral in ^31j, which indicates that the dimension of Q and Q is exactly 1/3. However, 
in the non-supersymmetric U{1) theories there typically is a lower bound for the conformal 
window. For example, in the extreme limit Np = we flnd the free Maxwell theory, which 
is not conformally invariant. We studied it on the of radius R in section 3 and found that 



-^Maxwell varics logarithmically with i?, eq. (3.25), indicating the lack of conformal invariance. 
One possible phenomenon for small Np is the chiral symmetry breaking in 3-dimensional 



QED coupled to massless fermions 16,17 . The numerical studies of lattice antiferromagnets 



51 suggest the QED theory with Nf = 8 Dirac fermions is a stable CFT, while the Nf 



theory is unstable to symmetry breaking towards a non-conformal ground state 10 . More 
generally, one of the signs of crossing the lower edge of the conformal window could be that 
the assumption of conformality leads to certain gauge invariant operators having scaling 
dimensions that violate the 3-dimensional unitarity bound A > 1/2. 

In 22,23 it was conjectured that F must be positive in a unitary CFT. Since as Np 
decreases so does F, it is possible that F may become negative for sufficiently small Np. 
This could serve as another criterion for theories outside the conformal window. It would be 
interesting to explore the different criteria above and to see if they are related. 
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A Vector spherical harmonics on 



In this appendix we review the vector spherical harmonics on S (see for example 52 ). 
The scalar spherical harmonics on transform in the [j, j] irrep of SU (2) x SU (2) for any 
half-integer j > 0. Denoting n = 2j + 1, the dimension of this representation is n^, and 
n > 1 is an integer. Let's denote these harmonics by Yngm, with < i < n and —i<m<i. 
Parameterizing by three angles (x, 0, 0), with the line element given by 

ds^ = dx^ + sin" X [de^ + sin' 6 d(j)^] , (A. 1 ) 

one can write down explicit formulas for the spherical harmonics: 

• £ 

sm v 

YnUx, 0, 0) = —^^ni{x)Yem{0, 0) , (A.2) 



where are the usual spherical harmonics on S"^ and 

^ , , (i^+^cos(nx) niT {i + n)\ 

The spherical harmonics satisfy 

{V^'V^ + n^~l)Ynem = 0. (A.4) 
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As mentioned in the main text, any one-form on can be written as a linear combination 
of forms that transform in S0{4) irreps with Jl = j^, ji = j/j + 1 and = ji + 1. One way 
to understand this fact is as follows. One can express any one-form as a linear combination 
of the right- invariant one-forms on S^. These right-invariant one-forms transform in the [1, 0] 
irrep. The coefficients of the right-invariant forms in the decomposition of an arbitrary form 
can in turn be expanded in terms of the usual spherical harmonics, which as mentioned above 
transform in the [j,j] irreps. The Hilbert space of square-integrable one-forms therefore 
decomposes as 



([1, 0] ® = ([j - 1, j] © © [j + 1, j]) 



(A.5) 



The sum runs over all j G N/2, but when j = 1/2 the first term in the paranthesis is absent, 
and when j = the first two terms are absent. Switching from j to n we see that the 
one-forms on transform in 



e 

71=2 I- 



1 n 



e 

n=2 I- 



n — 2 n 
'2 



e 

n=2 I- 



n n 



(A.6) 



< 



We will call the harmonics corresponding to the first sum SJ^*"^ (with < ^ < n and 
m < i, having total dimension ra^), and those corresponding to the second and third sums 
V^f™ and V^^™, respectively (with < i < n and —i<m<i, having total dimension n —1 
for either Yf^"^ or Y]l^). 

Explicit formulas are available. Defining the inner product on the space of square- 
integrable one-forms in the usual way. 



(A,: 



dx de del) sin^ X sin 9 A^(x, M^{x, 9, </>) 



(A.7) 



we normalize the harmonics so that 



) = (v: 
) = (§■ 



/\jnlm \m'£'m' 
R ) ^ R 



nim \Tn'£'m' 



) = (v: 



nim \Tn'i'm' 



) = 0. 



(A. 



As explained in section |3.2.2| the S"-^™- are gradients of the usual scalar harmonics: 



dYnlmiX^OA) 



(A.9) 
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and they are the only closed forms in the decomposition (3.15). The co-closed forms V 



jnlm 



and V^^'" can be recast into the symmetric and antisymmetric combinations 



V' 



nlm 



"^nlm _j_ "^nim 



jnlm 



•^nlm "^nlm 



V2 



(A. 10) 



which by virtue of (A. 8) are also orthonormal. Then 



sin^+^X^nKx) ,^ , 
^nH{i + l)a„, 



(A.ii: 



where *2 denotes the Hodge dual on S"^ with the standard line element. These expressions 
exhibit V"^*^ and W"^™ explicitly as co-closed one-forms. 



When we evaluate the sums in equations like (3.17) we need to know what happens to 



the one-forms close to the North pole at x = 0. A simple analysis of eqs. (A. 9) and (A.ll) 
gives: 



jnlm 



{x,o,4>)\ 



0{x''-') 
oix'') 



2l-2\ 



(A. 12) 



as X 0. The only harmonics that are non-zero at x = are therefore S"^™ and V"^™ with 
i = 1 and m = — 1, 0, 1. 
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